Testing the global star formation relation: An HCO + (3-2) 
mapping study of Red MSX sources in the Bolocam Galactic 

Plane Survey 

David E. Schenck 1,2 , Yancy L. Shirley 1,3 , Megan Reiter 1 , 

and 

Stephanie Juneau 1 
ABSTRACT 

We present an analysis of the relation between the star formation rate (SFR) 
and mass of dense gas in Galactic clumps and nearby galaxies. Using the bolomet- 
ric luminosity as a measure of SFR and the molecular line luminosity of HCO + 
(3-2) as a measure of dense gas mass, we find that the relation between SFR 
and Mdense is approximately linear. This is similar to published results derived 
using HCN (1-0) as a dense gas tracer. HCO + (3-2) and HCN (1-0) have similar 
conditions for excitation. Our work includes 16 Galactic clumps that are in both 
the Bolocam Galactic Plane Survey and the Red MSX Survey, 27 water maser 
sources from the literature, and the aforementioned HCN (1-0) data. Our results 
agree qualitatively with predictions of recent theoretical models which state that 
the nature of the relation should depend on how the critical density of the tracer 
compares with the mean density of the gas. 

Subject headings: Galaxies: star formation-Radio lines: general-Stars: formation 



Introduction 



All stars originate in clouds of cold, dense molecular gas (Evans 1999). The formation 
of low-mass stars has been studied in depth and the process by which it occurs is believed 
to be fairly well understood (Shu, Adams, & Lizano 1987, Andre et al. 2000). However, 
the formation of stars with masses greater than ~8 M is more difficult to examine because 
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the regions in which they form are considerably farther away (Wu et al. 2010). High- mass 
stars may not simply be a scaled-up version of low-mass stars, but instead may be formed 
in a fundamentally different way (e.g., competitive accretion, core coalescence vs. scaled-up 
accretion; see McKee & Ostriker 2007). A complete theory of star formation requires an 
understanding of the efficiency and rate of forming stars across the entire stellar mass range. 

In the last few years, a new nomenclature has been adopted to describe the structure 
of Galactic star-forming regions within molecular clouds. Protostars form from individual 
dense cores of gas. Most protostars form in clustered environments and collections of cores 
which may be physically connected or gravitationally bound within a larger clump. Clumps 
located in the Milky Way are near enough to be observed in detail using submillimeter 
single-dish telescopes with resolutions of tens of arcseconds (e.g., Reiter et al. 2011). The 
information gained from studies of these local star-forming clumps can be used to understand 
star formation in galaxies too far away to resolve individual star-forming regions. Ultimately, 
these studies may be used to constrain a universal star formation law. 

Determining the rate and efficiency of star formation has been a long standing problem. 
Schmidt (1959) first proposed that the rate of star formation would depend on the volume 
density of dense gas according to a power law. Building upon this assumption, it was 
shown globally for galaxies that the power law phrased in terms of surface density gave 
^sfr oc S^ ±tU5 (Kennicutt 1998). This relation, which was derived using CO to trace the 
molecular gas and HI to trace the atomic gas is called the Kennicutt-Schmidt relation (K-S 
relation). Even with the resolution of state-of-the-art interferometers, the K-S relation has 
only been probed down to scales of ~ 1 kpc, much larger than the sizes of giant molecular 
clouds (tens of pc) in nearby galaxies (e.g., Bigiel et al. 2011). The K-S relation has 
been questioned because CO observations probe low density molecular emission and not 
the dense molecular gas more directly associated with the clumps and cores from which 
protostars form (e.g. Gao & Solomon 2004a). Therefore, subsequent studies have tested 
the power law relation by measuring the molecular line luminosity of tracers that specifically 
trace the dense gas mass. Infrared luminosity is used as a measure of the star formation rate; 
collapsing protostars emit strongly in the IR. Using the line luminosity of the transition HCN 
(1-0) and infrared luminosity of a set of 65 galaxies, including luminous and ultraluminous 
infrared galaxies (LIRGs and ULIRGs), Gao & Solomon (2004a,b) found a power law index 
of 1.00 ± 0.05. They interpreted this result as the SFR being proportional to dense gas mass 
(Gao & Solomon 2004a,b). This HCN (1-0) survey was extended to the scale of Galactic 
dense clumps with similar results; the power law index was very near one over a range of 
7 orders of magnitude in L IR (Wu et al. 2005, Wu et al. 2010). This was interpreted as 
evidence that there is a fundamental unit of star formation associated with massive clumps 
with luminosities L& ^ > 1O 4 ' 5 L and the star formation efficiency within the dense clumps 
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is constant. 

Since the initial work of Gao & Solomon and Wu et al., a wider range of critical densities 
has been probed by different molecular tracers (e.g., Bussmann et al. 2008, Juneau et 
al. 2009). Theoretical radiative transfer models predict the behavior of the SFR and line 
luminosity as a function of gas density on galactic scales (Krumholz and Thompson 2007, 
Narayanan et al. 2008). The behavior of the line luminosity is dependent on the critical 
density of the tracer in these models; if the critical density exceeds the mean density, the line 
luminosity rises superlinearly with density. However, for critical densities below the mean 
density, it rises linearly. The SFR is assumed to rise superlinearly with density in these 
models 0, similar to the original Schmidt relation (1959) but different from the conclusion 
of Gao & Solomon. Combining the two conditions gives linear correlations for high critical 
density tracers and superlinear correlations for low critical density tracers. This general 
trend is observed in galaxies using CO (1-0), HCN (1-0) and (3-2), and HCO + (1-0) and 
(3-2) (Juneau et al. 2009, also see Baan et al. 2008, Iono et al. 2009). 

Gao and Solomon (2004) and Wu et al. (2010) interpret their results as evidence that 
the SFR to dense gas mass ratio is constant. However, the extragalactic community and 
many theoretical models continue to assume a superlinear Schmidt relation. In this paper, 
we test the observed Gao & Solomon and Wu et al. correlation by mapping a sample of 16 
massive clumps in HCO + (3-2) selected from the Bolocam Galactic Plane Survey (Aguirre 
et al. 2011). This sample is combined with 27 massive clumps mapped in HCO + (3-2) that 
are associated with water masers from Reiter et al. (2011) and 14 (U)LIRGs, from Juneau 
et al. (2009) that were adapted from Gracia-Carpio et al. (2008) using updated distances. 
We also re-analyze the 50 dense clumps from Wu et al. (2010) and 42 galaxies from Gao 
and Solomon (2004) to verify the published HCN correlation. 

We choose HCO + (3-2) because of its ability to trace nearly identical excitation con- 
ditions to HCN (1-0). The frequency and dipole moment of HCO + (3-2) are a factor of 
~3 and ~1.5 times higher, respectively, than the corresponding values of HCN (1-0). The 
Einstein A coefficient is proportional to z/ 3 /z 2 , making the Einstein A coefficient of HCO + 
(3-2) ~ 60 times higher. Most of this disparity is compensated by HCO + inducing dipoles 
in the molecules with which it interacts, increasing the effective cross section for collision 
and the corresponding collision rate 7^. HCN is not an ion, so it does not induce a dipole 
in nearby molecules and it has a comparatively lower rate of de-excitation. At 20 K, HCO + 
(3-2) and HCN (1-0) have critical densities of 3.90 x 10 6 and 1.25 x 10 6 cm -3 , respectively, 
a factor of ~3.1 difference. Since the excitation of a molecular line varies gradually over a 
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range of densities, claiming the line traces gas at the critical density is an oversimplification 
(Evans 1999). A more practical quantity is the effective density for excitation, n e ff, which 
is defined as the density at which a transition will have a radiation temperature of 1 K 
assuming log(N/Av) = 13.5 cm _2 /km/s and an average kinetic temperature of 20 K (Evans 
1999, Reiter et al. 2011). The two transitions also have even more similar effective densities; 
HCO+ (3-2) and HCN (1-0) have effective densities of 2.33 x 10 4 and 1.50 x 10 4 cm" 3 . These 
only differ by a factor of ~1.5. Densities can range over several orders of magnitude within 
molecular clouds, so such a small difference means that the transitions essentially trace the 
same molecular gas in the absence of strong chemical differentiation. 

In §2, we summarize the observations of our sample of 16 massive cores as well as the 
full (combined) sample used in this paper. In §3, we discuss the methods used to derive 
physical quantities from the maps. We present our results and compare them to both other 
observations and theoretical predictions in §4 as well as critically discuss several caveats with 
the interpretation of the observed correlation. 

2. Observations 

Sixteen sources were mapped using the 10-meter Heinrich Hertz Telescope (HHT) on 
Mount Graham, Arizona. Observations were made on April 10 and 11 as well as May 4 
through May 9, 2010, excluding May 5. The HHT was equipped with the 1mm ALMA 
prototype sideband-separating receiver. The observations were made with a 6.0 GHz IF in 
dual polarization 4 IF mode with the vertical and horizontal polarizations each split into 
upper and lower sidebands. Two transitions were observed simultaneously; HCO + J = 3 — > 2 
(267.558 GHz) was centered in the lower sideband while N 2 H+ J = 3 ->• 2 (279.512 GHz) was 
centered in the upper sideband. The N 2 H + maps are not included in this analysis because of 
the weakness of the detections. The observed line properties are listed in Tables 2 and 3. The 
backend filter-banks were split into 256 channels with 250 kHz resolution. At the frequencies 
used, the velocity resolution was ~0.28 km s -1 . The maps have an angular resolution of 
27.2". 

Maps were made using the On-The-Fly (OTF) imaging technique (Mangum et al. 2007). 
Scan rows were separated by about 10", about one-third of the beam size. Each source was 
scanned at least once in the RA and DEC directions and many were scanned again in both. 
Additional pairs of scans were offset by about 5" from the first pair to better sample the 
source. The pointing accuracy of the HHT is typically 5" RMS. 

The maps were reduced using the GILDAS CLASS software. First, a linear baseline 
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was removed from each map. The vertical and horizontal polarizations were scaled sep- 
arately using their respective main beam efficiencies and then combined (?7#=0.70±0.04, 
77y=0.81±0.04). The maps were convolved using a Gaussian-tapered Bessel function of the 
form ^Mexp[-(D 2 ] where a = 1.55(^f ) and b = 2.52(^f) (Mangum et al. 2007, Reiter 
et al. 2011) to preserve spatial resolution. 

2.1. Source Selection 

The sixteen sources mapped for this project were selected from the Bolocam Galactic 
Plane Survey (BGPS). A total of 1402 BGPS sources spectroscopically detected by Schling- 
man et al. (2011) were examined by eye using GLIMPSE (8 /iin) and MIPSGAL (24 /zm) 
images. The MIPSGAL images were used to determine which sources were associated with 
a 24 /iin point source. The point sources indicate the presence of either an embedded proto- 
star or an evolved star. Each BGPS source within 15" of the point source was considered to 
contain a candidate embedded object. The 247 objects associated with 24 fj,m point sources 
were then compared to the Red MSX Source (RMS) survey by calculating the offset between 
their positions. The RMS survey includes over one thousand massive young stellar objects 
(MYSOs) that are radio-quiet, mid-infrared point sources (Mottram et al. 2010). Comparing 
the coordinates of the objects in each catalogue, 79 from MSX were within 15" of a BGPS 
source observed by Schlingman et al. The Schlingman et al. subsample of BGPS sources 
contains the highest flux 1.1 mm sources in the BGPS. Candidate RMS sources in the BGPS 
were observed generally in order of brightest to dimmest HCO + (3-2) intensity. Not every 
object common to the two surveys was observed due to observational constraints, but every 
joint object with Ihco+ (3-2) > 10.5 K km s _1 was mapped. 

In addition to our observations, we include HCO + (3-2) maps from the survey of Reiter 
et al. (2011). Twenty-seven massive clumps associated with H 2 masers were taken from 
this paper. The luminosity and distance of the BGPS-RMS sources and water maser sources 
are listed in Table 1. The bolometric luminosity is given for the BGPS-RMS sources while 
the infrared luminosity is given for the water maser sources. Also, fourteen luminous and 
ultraluminous infrared galaxies (LIRGs and ULIRGs) originally from Gracia-Carpio et al. 
(2008) represent the currently known extragalactic observations of HCO + (3-2). We analyze 
the correlation between bolometric luminosity and molecular line luminosity for this com- 
bined sample using a Bayesian linear regression technique (Kelly 2007). We also re-analyze 
the Wu et al. (2010) correlation of 50 massive, dense clumps and 42 galaxies that were 
mapped in HCN (1-0) (Wu et al. 2010, Gao and Solomon 2004a,b) using the same Bayesian 
linear regression technique. This was done in order to characterize the uncertainty in the 
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slope of the correlation as well as make the correlations of the two tracers more directly 
comparable. Galaxies which had upper or lower limits on the molecular line luminosity were 
not used in our analysis. 



3. Clump Physical properties 

The maps in integrated intensity of HCO + (3-2) emission are shown in Figure 1. Inte- 
grated intensities are given by, 

I(T mi ) = jT mb dv±VSv Bne Sv ckan a^ (1) 

where 8vu ne is the velocity extent of the line, 5v c h an is the channel width of the spectrometer 
(0.28 km s" 1 ), and OT mb is the uncertainty in the main-beam temperature. The quantity 
Svu ne was derived from setting boundaries on each side of the line so that the entire line was 
enclosed. As a result, the intensity uncertainty is a conservative estimate. The line with 
the highest integrated intensity in each map was used in calculating other quantities. The 
peak integrated intensities are usually near the center of the map, but that is not always 
true, especially when there are two objects in a single map. The total sample of galactic 
cores mapped in HCO + (3-2) had a mean intensity of 52.6 ± 50.2 K km -1 and a median of 
36.5. The standard deviation of the sample is relatively large because the sources cover a 
wide range of intensity from 3.2 to 195.9 K km s -1 . The mean and median for the subset 
of BGPS/RMS selected cores are 13.6 ± 9.5 K km s^ 1 and 11.7 K km s _1 . The mean and 
median for the subset of water maser selected sources are 75.6 ± 50.3 K km s^ 1 and 61.5 K 
km s _1 (Reiter et al. 2011). The BGPS selected subset extends the Reiter et al. sample of 
Galactic sources to lower integrated intensities in HCO + (3-2). 

The angular diameter was chosen to be the area within the the half-peak intensity 
contour, A (6 s *b = y/AApK). The uncertainty in the half-peak diameter was determined 
by finding the diameters associated with the half-peak contour ± c/(T m6 )- The beam of 
the telescope (27.2") contributed to the width of the maps. This effect was removed by 
calculating the deconvolved sizes using 



ddec ~ y #L - d leam ( 2 ) 

The distance and angular diameter were used to calculate a physical radius from 

R = ^ (3) 

where R is the radius in the units of the distance, D, and 9d ec is the deconvolved angular 
diameter in radians. The total sample has a mean size 0.31 ± 0.18 pc and a median of 0.28 
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pc. The means of the BGPS/RMS and water maser samples are 0.35 ± 0.19 and 0.29 ± 0.16 
pc and medians 0.34 and 0.24 pc, respectively. The BGPS/RMS mean is skewed upward by 
a large source (019.884-053) which is also the brightest in that sample. A histogram of the 
sizes is presented in Figure 2 along with histograms of other physical properties. 

Distances were obtained by using radial velocity measurements to get kinematic dis- 
tances (Schlingman et al. 2011). Kinematic distances are ambiguous in that two different 
distances correspond to the same radial velocity in the first quadrant of the Galaxy. Four 
sources were assigned to the near distance because they appear to be associated with an 
infrared dark cloud (IRDC) as determined by visual inspection of MIPSGAL images. The 
remaining distances were determined using information supplied by the RMS online catalog 
(Urquhart et al. 2008) which not only gives kinematic distances, but also a "complex dis- 
tance." The complex distance was for any MYSOs belonging to a complex with a previously 
determined distance. If one of the two kinematic distances was within about 1 kpc of the 
quoted complex distance, that kinematic distance was assigned. The water maser clumps 
from Reiter et al. (2011) used distances that were previously determined (see Shirley et al. 
2003). 

The uncertainty in the distances is needed to propagate the uncertainty in the molecular 
line luminosity. The uncertainty was determined using distance probability curves produced 
by Ellsworth-Bowers et al. (in prep.). This method employs radial velocity measurements, 
proximity to IRDCs, H I self- absorption, and other factors to calculate a joint probability 
distribution (P(D)). The distance probability curves are typically bi-modal (centered on the 
near and far kinematic distances). Each mode of the probability curves were fit by Gaussians, 
one for the near distance and one for the far distance, to determine the spread of the curves. 
Probability curves were available for all of the BGPS/RMS crossover sources, but only 8 of 
the 27 water maser sources (Ellsworth-Bowers, private communication 2010). The remaining 
19 sources without distance probability curves were assigned the average uncertainty of those 
8. There was no trend between distance and distance uncertainty. 

The bolometric luminosities were calculated using the observed spectral energy distri- 
butions (SED) of each source. SEDs were constructed uniformly from fluxes taken from 
the BGPS (1.1 mm), MSX (8, 12, 14, and 21 fim), SCUBA (850 /im), MIPSGAL (70 /im), 
IRAS (12, 25, 60, and 100 /im), and AKARI (65, 90, 140, and 160 /im). Bolometric lumi- 
nosities were calculated by integrating over the observed fluxes. Luminosities were already 
available for the water maser sources, but they were re-calculated because new fluxes had 
become available since the published SEDs in Mueller et al. (2002). See Table 4 for the 
complete SEDs of the sources. The logarithm of the mean luminosity (in L ) of the sample 
is 5.09 ± 0.87 with a median of 4.52. The mean bolometric luminosity of the water-maser 
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sources exceeds that of the BGPS/RMS sources; the former is 5.23 ± 0.76 and the latter 
4.64 ±0.66. 

L', or the molecular line luminosity, is a commonly used quantity in extragalactic studies 
as an analogue to the mass of molecular gas. Formally, molecular line luminosity, or source- 
integrated surface brightness as it is called in Mangum et al. 2008, is defined as 

U = 23.5047^(1 + z)- 3 D 2 L I(T mb ) (4) 

where 9 s *b is the convolved angular diameter in arcseconds, z is the redshift, is the 
luminosity distance in Mpc, and I{T m ^) is the integrated intensity in K km s -1 . Two sources 
with the same main beam brightness temperature and spatial extent will have the same L'. 
The mean of the line luminosity for the BGPS/RMS and water-maser samples were each 33.5 
K km s -1 . Histograms of the bolometric and line luminosities are shown in Figure 3. In §4.1 
we construct the correlation between bolometric luminosity and molecular line luminosity 
with distance uncertainties propagated for the Galactic sources. 

Finally, we calculated the mass in the Galactic clumps using two different methods. 
Virial masses were calculated using 

M mr = 5R(AV)2 « 209 ^ A ^ 2 (5) 
8aid2Gln(2) a±a2 

1-P/3 fR . 

ai = Y~^l (6) 

where R is the physical radius in pc, Av is the FWHM linewidth in km s -1 , and a\ and a<i 
are factors that account for a power law density distribution (n oc r~ p ) and non-spherical 
shape, respectively (Bertoldi and McKee 1992). The median value of p (1.75), the exponent 
of the density power law, was adopted from Mueller et al. (2002), making ai=1.39 for the 
BGPS/RMS crossover sources. The values of a\ were determined in Mueller et al. (2002) 
for the water-maser sources. The correction factor a 2 is negligible for aspect ratios below 
2, so it was set to 1 for all calculations. The linewidth is a measure of the thermal and 
turbulent support against gravitational collapse, but can be broadened by optically thick 
emission, leading to overestimates of the virial mass (Shirley et al. 2008, Reiter et al. 2011). 
The logarithm of the mean virial mass (in M Q ) is 3.30 ± 0.41 and a median of 3.13. For the 
BGPS/RMS and water-maser samples, the mean is 3.24 ±0.41 and 3.33 ±0.41, respectively. 

Masses were also estimated from the 1.1 mm continuum emission by assuming a single 
dust temperature and common dust opacity for all sources. The isothermal mass is 
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where is the total integrated flux at 1.1 mm and k v is the opacity of that wavelength. The 
value /ti.i mm =0.0114 cm 2 g" 1 was adopted (Ossenkopf & Henning 1994, Enoch et al. 2006, 
Battersby et al. 2010), consistent with a dust-to-mass ratio of 100. Assuming a temperature 
of 15 K, the temperature used in Battersby et al. (2010), the equation for M iso becomes 

M, so = 14.32(e^ - 1) (A) (^)V (8) 

The logarithm of the mean isothermal mass (in M Q ) is 3.37 ± 0.64. This is higher than 
the mean for the virial mass. The water-maser sample also has a higher mean isothermal 
mass (3.54 ±0.51), but the mean isothermal mass is lower than the mean virial mass for the 
BGPS/RMS sources (3.09 ± 0.87). 



4. Results and discussion 

A Bayesian linear regression routine from Kelly (2007) was used to test the correlation 
between L' and L bo i as well as L' and mass. The program uses a Markov chain Monte 
Carlo (MCMC) method to take random draws of the slope, intercept, and intrinsic scatter 
by perturbing the previous set of line parameters. These new sets of parameters are either 
accepted or rejected according to the Metropolis-Hastings algorithm (Metropolis et al. 1953, 
Hastings 1970). The program tests 10,000 sets of parameters and saves the distribution of 
accepted values. The spread in the distribution of the slopes is entirely attributable to the 
uncertainties in the data; the uncertainty from the MCMC sampler is negligible (Brandon 
Kelly, private communication 2011). For a more detailed description of the program, see 
Kelly (2007). 

The results of the Bayesian linear regression fits are presented in Table 5. The lines 
take the form 

log(y) = b x log(a;) + a. (9) 

The ordinate and abscissa variables are listed in Table 5 in the order y-x. The quantity Oi nt 
is the intrinsic scatter of the fit. 



4.1. M-L/ correlation 

First we test the correlation between mass and line luminosity. The masses of Galactic 
clumps were found using both isothermal mass (determined from dust continuum emission) 
and virial mass (determined from HCO + emission). The Mj SO -L' relation (Figure 3(a)) had 
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a sublinear slope (0.78 ± 0.10) and the M vir -L' relation (Figure 3(b)) had a linear slope 
(1.06 ± 0.18). Using M iso produced a slightly tighter correlation; the correlation coefficient 
for M iso was 0.85 compared to only 0.70 for M vir . The same trends were observed in multiple 
molecular species in Reiter et al. (2011) and Wu et al. (2010) using the water maser sources. 

The virial mass depends on the linewidth, Av. The line can be artificially broadened if 
the emission is optically thick which is true for the transitions considered (see Phillips et al. 
1979). This effect will cause the virial mass to be an overestimate. M iso , on the other hand, 
does not depend on the line parameters, unlike M vir and L' . Thus, we prefer M iso as it is 
an independent estimate of the clump mass. Nevertheless, our tight linear correlation with 
both mass estimates indicates that L' is indeed a good tracer of the dense gas in clumps. 



4.2. L'-Lboi correlation 

Correlations were found for several subsets of the data by using L' as the independent 
variable and as the dependent variable. These are presented in Figure 3(c) (HCO + 
(3-2)) and Figure 3(d) (HCN (1-0)). The slope of the galactic clumps and combined data 
sets observed in HCO + (3-2) both came out to nearly one (1.10 ± 0.12 and 1.04 ± 0.02 
respectively). The histograms in the bottom right corner of each panel in Figure 3 show the 
distribution of slopes for the galactic and extragalactic data. The slope for all clumps and 
galaxies together has an uncertainty nearly a factor of 10 lower than the clumps alone. This 
is the result of the data appearing in two widely separated groups in L^-L' space so that the 
line was basically fit to two clumps of points. The slope of the correlation for Galactic clumps 
agrees very well with the Wu et al. (2010) correlation for HCN (1-0) (also see below). While 
the best-fit line for the Galactic clumps extends to be slightly higher than the extragalactic 
objects, the uncertainty on the slope makes it such that the extragalactic points are within 
la of the fit. The extrapolation of the best-fit of the extragalactic points misses the Galactic 
points more dramatically, although the large uncertainty on the slope still cannot rule out 
a linear slope at the more than la level. This result is different than was found for HCN 
(1-0) by Wu et al. (2005), however, we must caution that the published extragalactic HCO + 
(3-2) detections are small in number and limited to > 10 11 L Q . Global extragalactic 
observations of L' in HCO + (3-2) are needed, especially for galaxies with Lj, i < 10 11 L to 
better compare with the HCN (1-0) results. 

Linear fits for the HCN (1-0) data set were presented in Gao and Solomon (2004a,b) 
and Wu et al. (2010), both of which found a linear relation between Lir and L'. We test 
their results by re-analyzing their data with the Bayesian routine. The line luminosity of one 
source, DR21S, was misprinted in Wu et al. (2010) (it is 2.29 rather than 0.29; Wu, private 
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communication 2011). This did not affect their results as their correlations were found using 
the correct value. The fits displayed in the top right plot of Figure 2 do not include any 
points below L IR =10 3 L Q . This cutoff was tested because that is roughly the lowest L bo i in 
the HCO + data set. The limit makes the HCN (1-0) and HCO + (3-2) correlations span the 
same range in luminosity. The slopes for the galactic, extragalactic, and combined sets are 
1.17 ± 0.11, 1.07 ± 0.06, and 0.99 ± 0.01 respectively. These slopes are all within 1.5 a of 
unity. The agreement between the extrapolated galactic and extragalactic best-fit lines is 
better than for HCO + (3-2). The extragalactic correlation from Gao & Solomon (2004a,b) 
spans 1.5 more orders of magnitude in Lw than the observed HCO + (3-2) in galaxies and 
likely provides a more robust sample from which to determine the fit. Our re- analysis of the 
Gao & Solomon and Wu et al. data confirms their results. 

The amount of dense gas in the galactic clumps needs to be high enough to fully sample 
the IMF. The IMF describes the number of stars per mass bin and, generally, lower mass 
stars are more common. A cloud with a relatively small mass will not form as many stars 
which makes the probability of seeing a massive star lower. When we say " sample the IMF" , 
we mean the probability of seeing stars across the full range of the IMF is high. Based 
on the assumption that larger clouds with more stars will be more luminous, Wu et al. 
(2010) instituted a minimum Ljr of 10 4 ' 5 L Q . The linear regression was performed with this 
condition imposed in order to directly compare to their work. The slopes were consistently 
lower; the galactic and combined samples had slopes 0.88 ± 0.15 and 0.94 ± 0.02, respectively. 
These slopes may be smaller because of the way the cutoff was imposed. Simply removing all 
points below a certain luminosity will allow more scatter in the horizontal direction compared 
to the vertical which will flatten the line. 

In this work, we choose to define the cutoff as a minimum line luminosity rather than 
infrared or bolometric luminosity. The motivation behind the cutoff is to define a minimum 
mass which we have constrained from the observed correlations between line luminosity and 
mass. Krumholz and Thompson (2007) claim that the cloud mass should exceed ~1000 M 
to accomplish a decent sampling. Using the correlation between L' and Mj S0 , setting the 
mass to 1000 M Q gives a line luminosity of 0.85 K km s -1 . Based on the Galactic L^-L' 
correlations with no cutoff, this corresponds to bolometric (or infrared) luminosities of 10 4 ' 2 
and 10 3 ' 4 for HCO + and HCN, respectively. The region containing points that are below 
both cutoffs are indicated by dotted lines in the bottom left corner of Figure 3(c) and 3(d). 
In both cases, most of the points in these regions lie below the correlation, indicating they 
may not be sampling the IMF well enough. There are not many points in these regions and 
the downturn observed in HCN (1-0) is not apparent in HCO + (3-2), so it is not possible 
to definitively conclude where the cutoff should be, but there is no evidence against those 
derived assuming 1000 M as the minimum mass necessary. 
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4.3. Comparison to theory 

The physical basis for the observed L'-L^ correlations have been studied using numerical 
simulations. Krumholz and Thompson (2007) modeled isothermal and homogeneous giant 
molecular clouds (GMCs) using a radiative transfer code and considering star formation 
regulated by turbulence. Narayanan et al. (2008) used 3D non-LTE radiative transfer with 
hydrodynamic simulations of isolated galaxies and galaxy mergers. Despite considering two 
different classes of objects, the same explanation was reached by both. The origin of the 
SFR-L' relation comes from the Schmidt relationship and the density dependence of L'. The 
Schmidt relation says that the star formation rate is proportional to density to some power, 
or 

SFR oc p N (10) 

where N is the Schmidt index. The value of N is 1.5 with one factor of p coming from 
the mass available for stars and a factor of p ' 5 coming from the density dependence of 
the dynamical time for collapse (Madore 1977, Elmegreen 1994). This assumes that star 
formation occurs at a constant rate per dynamical time. To predict the relation between 
SFR and line luminosity, 

SFR(xL' a , (11) 
the dependence of the line luminosity on density must be known. This is expressed as 

l! oc p p (12) 

so that the index of the SFR-L' correlation is oc—j. 

Both simulations determined that density dependence of L' fell into two regimes, one in 
which the critical density is much larger than the mean density of the gas and one in which 
it is much smaller. When n crit » n, only the high density tail of the density distribution 
is above the critical density so emission predominantly arises from the densest regions. The 
emission does not trace the density of the gas, but instead will trace the same density in 
every subcritical object. As density increases (with n crit remaining much larger than n), the 
line emission rises superlinearly so that /3 > 1 and a < N. The SFR-L' relation will have an 
index below the Schmidt index. Both HCO + (3-2) and HCN (1-0) fall in this regime. In 
the other case where n crit « n, nearly all of the gas is thermalized and emitting. The line 
luminosity faithfully traces the mass and rises linearly with density. The index of the SFR-L' 
relation will approximately equal the Schmidt index. CO (1-0) produces an index near the 
Schmidt index, which is consistent with this formulation as CO traces lower densities than 
HCO+ and HCN. 

Some combination of gas conditions and molecular transitions will be intermediate be- 
tween the two cases. As density increases, the relation passes from the subcritical to the 
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supercritical case, meaning the SFR-L' relation should transition from linear to superlinear. 
Evidence of this upturn have been found in observations, though it is only tentative due to 
a caveat to be discussed in the next subsection (Gao et al. 2007). 

A power law index of 1.0 was predicted for HCN (1-0), in agreement with observations. 
The critical density of HCO + (3-2) is about a factor of 3 higher than that of HCN (1-0). 
Gas tracers have critical densities that span several orders of magnitude, so the difference 
in n crit for the two tracers considered here is relatively small and they can be considered 
to trace the same conditions. Based on the theoretical predictions, the index for HCO + 
(3-2) should also be near 1.0. The results presented in this paper are consistent with the 
theoretical predictions of Krumholz & Thompson and Narayanan et al. The interpretation 
that a linear relation means a constant ratio between the SFR and dense gas mass cannot 
be distinguished from the theoretical predictions based on our work. 



4.4. Caveats 

The Lboi-L' relation is intended to represent a relation between the SFR and gas mass. 
However, the conversions are complicated by a wide range of assumptions which must be 
explored. 

Bolometric (or infrared) luminosity is the most popular measure of the star formation 
rate (SFR) as it is easy to determine. On scales from 1 kpc to galactic scales, the IMF of stars 
is well sampled and the conversion from L boi to SFR is linear (see Kennicutt 1998). However, 
this conversion for galaxies is potentially confounded by the presence of active galactic nuclei 
(AGN) which contribute their own IR emission by heating dust grains (Andreani et al. 2003, 
Carilli et al. 2005, Li et al. 2008). The extra observed IR flux can lead to an overestimate 
of the SFR. 

The effect of AGN on the L bor L/ correlation is expected to be negligible. For a sample 
of 68 ULIRGs, about 85 % of the infrared luminosity was found to be powered by star 
formation with only 15 % attributable to AGN (Nardini et al. 2008). While AGN occur 
more frequently in more IR-luminous galaxies (e.g., Veilleux et al. 1995), their fractional 
contribution to the total bolometric luminosity mostly remain small. The net effect on the 
Lboi-L' relation was tested in Juneau et al. (2009) by comparing the slopes obtained when 
using either the total infrared luminosity (L IR ) or the far-infrared luminosity (L FIR ) as a 
proxy for L^o/. The far-infrared luminosity is expected to be dominated by star formation 
even in the presence of AGN-heated dust emission, which peak in the mid-infrared. In 
contrast, the total infrared emission is more subject to AGN contamination. However, the 
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Lb i-L' relation is unchanged when substituting for Lpm, suggesting that possible AGN 
contributions do not affect the relation (at least within the uncertainties). 

While Galactic clumps will not be affected by AGN, their L hoi also is not a perfect 
measure of the SFR. The most massive stars in a system younger than ~3 Myr will not 
have had time to evolve and die in a supernova; therefore, not enough time has elapsed to 
achieve equilibrium between stellar birth and death. In this situation, the light emitted will 
not trace the star formation rate. A simulation run by Krumholz and Tan (2007) found 
that the luminosity per unit star formation is insensitive to the age for populations older 
than a few Myr (which is a timescale that is easily averaged over 1 kpc within a galaxy but 
would not necessarily apply on the scale of an individual Galactic clump). The luminosity 
of younger populations traced mass, but not the star formation rate. Dividing the mass by 
an independent measure of the age to estimate the star formation rate does not solve the 
problem for systems younger than ~1 Myr because not all of the stars in such a system 
have finished collapsing. These pre-main sequence stars are powered by accretion and will 
release more energy than a main sequence star of the same mass. Since the contribution of 
these collapsing objects is unknown, it is impossible to determine how much mass is in the 
system (Krumholz and Tan 2007). By resolving the infrared population in Galactic clumps 
(e.g., with James Webb Space Telescope or Large Binocular Telescope Interferometer) and 
directly comparing the mass distribution within the clumps (e.g. with ALMA), we may be 
able to better constrain the actual SFR in Galactic clumps. 

Chemical variations potentially make the utility of HCN and HCO + questionable as 
uniform dense gas tracers in different environments. The strong X-ray emission from an AGN 
can increase ionization and therefore the abundance of HCN in relation to other molecules 
(Lintott & Viti 2006). This is the explanation given for the trend observed in the HCN(1- 
0)/HCO + (1-0) luminosity ratio which increases with increasing far-infrared luminosity for 
Lf/r > 10 11 L by Gracia-Carpio et al. (2008). However, another interpretation is that more 
luminous galaxies have a higher dense gas fraction and the effects on HCN chemistry are 
negligible. Using 5 transitions from 3 molecules (HCN, HCO + , CO) to form 10 line luminosity 
ratios (with the higher critical density tracer in the numerator), Juneau et al. (2009) showed 
that the ratios consistently increased with L IR except for the HCO + (3-2)/HCN(l-0) ratio 
which combines the two tracers considered here. This ratio, which has tracers of nearly the 
same critical density, was approximately constant. Using the same simulation as Narayanan 
et al. (2008), the line ratio behavior was replicated without including chemical abundance 
changes. They concluded the dense gas fraction increased with luminosity and that inclusion 
of chemical effects in the dense gas were not necessary. 

Molecular gas tracers are supposed to indicate the presence of molecular gas through 
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collisions, emitting more energy when they undergo more collisions. Ideally, all emission 
would be a result of interactions with the gas, but the emission can be affected by the 
background radiation. The excitation of HCN lines may be enhanced by IR pumping of a 14 
fim vibrational transition (Gracia-Carpio et al. 2006). Similarly, HCO + excitation can be 
affected by a 12 vibrational transition. Any emission associated with radiative pumping 
will lead overestimates of the gas mass inferred from the line luminosity. 

The complications discussed above present real obstacles to interpreting the L^-L' 
relation. None of them were included in the GMC or galaxy radiative transfer modeling 
because of the high degree of complexity that each would require. It is argued in both 
Krumholz and Thompson (2007) and Narayanan et al. (2008) that their ability to replicate 
the relations for different tracers without accounting for these details shows that their effects 
may be negligible. Despite these caveats, the results from this study show similar trends for 
the L bol -L' relation of HCO+ (3-2) as was observed by Wu et al. (2005, 2010) for HCN (1-0). 



5. Summary 

We conclude that line luminosity in HCO + (3-2) is tracing mass as evidenced by the fits 
of L' with isothermal mass and virial mass. M vir rises roughly linearly with L' while M iso 
also increases with L', but with a sub-linear slope. 

The slope of the log(L 6oi )-log(L') relation for HCO+ (3-2) is 1.04 ± 0.02 (1.09 ± 0.12 for 
Galactic clumps and 0.81 ± 0.21 for galaxies). This would imply a linear relation between 
the star formation rate and the mass of dense gas, though this interpretation is subject to 
several caveats. A similar result was seen by Gao and Solomon (2004) and Wu et al. (2010) 
using HCN (1-0), a tracer with similar conditions for excitation as HCO + (3-2). Those 
results were tested using the Bayesian linear regression routine. Considering all sources with 
L IR exceeding 10 3 L , the slope is 0.99 ± 0.01 (1.17 ± 0.11 for Galactic clumps and 1.07 ± 
0.06 for galaxies). 

The results of our work are generally consistent with theoretical predictions from radia- 
tive transfer models. Two different models, one considering GMCs (Krumholz and Thompson 
2007) and one considering isolated galaxies and galaxy mergers (Narayanan et al. 2008), pre- 
dicted that the slope of the log(L(, ;)-log(L') relation should depend on the critical density 
of the tracer and the mean density of the gas. If the critical density is above the mean 
gas density, then only high density peaks are thermalized and emit. In this case, the line 
luminosity rises super- linearly with density and the log(Lb ;)-log(L') relation is linear. If the 
critical density of the tracer is below the mean gas density, then nearly all the gas is traced 
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and the line luminosity faithfully traces the gas, rising linearly with density. This leads to 
a super-linear Lfe ;-L' relation. Both models predict a roughly linear relation for HCN (1-0) 
and since HCO + (3-2) is excited in similar conditions, it should behave similarly. 

The conversion of L bo i to SFR has many caveats. It can be confounded by the presence of 
AGN boosting the luminosity by heating dust grains or radiatively pumping the gas tracers 
to cause excess line emission. Uncertainty in the age of the stellar populations and in the 
ability to sample the IMF plague Galactic clumps. Ultimately, high resolution mid-infrared 
and submillimeter observations are needed to resolve embedded infrared populations and 
dense gas cores in Galactic clumps to reliably determine the SFR. All of these caveats must 
be considered and characterized before the underlying star formation relation can be fully 
understood. 

Our study suffered from a small number of galaxies mapped in HCO + (3-2) as only 14 
were available (Gracia-Carpio et al. 2008). A more complete survey must be conducted to 
extend the sparsely populated moderate luminosity region of the L^-L' plot toward < 
10 11 L . A followup extragalatic survey should include galaxies that span a wider range in 
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3.7 + 0.3 




40.50+2.54 


18 56 10 


7 53 14 


4.29 + 0.05 


2.1 + 0.4 




35.20-0.74 


18 58 13 


1 40 36 


4.76 + 0.05 


3.3 + 0.4 




59.78+0.06 


19 43 12 


23 43 54 


4.02 + 0.03 


2.2 + 0.6 




ONI 


20 10 9 


31 31 37 


5.11 + 0.02 


6.0 + 0.4 
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Table 1 — Continued 



Source 


a (J2000) 


<5 (J2000) 


log(L boi or L IR ) 


Distance Distance 1 




[h m s] 


[° ' "] 


[Lo] 


[kpc] Comment 


ON2S 


20 21 41 


37 25 29 


5.59 ±0.02 


5.5 ±0.4 


W75N 


20 38 37 


42 37 37 


5.29 ±0.03 


3.0 ±0.4 


W750H 


20 39 1 


42 22 50 


4.40 ±0.02 


3.0 ±0.4 


S140 


22 19 18 


63 18 47 


4.34 ±0.02 


0.9 ±0.4 


CepA 


22 56 18 


62 1 46 


4.52 ±0.02 


0.73 ±0.42 


NGC7538-IRS9 


23 14 2 


61 27 20 


4.59 ±0.05 


2.8 ±0.7 


S157 


23 16 4 


60 1 41 


4.31 ±0.03 


2.5 ±0.7 


Mean 






5.09 


4.11 


Standard Deviation 






0.87 


2.67 


Median 






4.52 


3.5 



1 IRDC: The source is associated with an infrared dark cloud. The near kinematic distance was assigned. 
Complex: The source is a part of a complex with a known distance. If one of the kinematic distances was 
similar to the complex distance, that kinematic distance was assigned (Mottram et al. 2010). 
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Table 2. HCO+ Properties 



Source 


T 
1 


VLSR 


v s*b 


17 
XX 


log(L ) 




[K km s i ] 


[km s 1 \ 


17/1 

[ ] 


r i 
[PC] 


[K km s 1 pc \ 


019.884-053 


46.12 + 2.47 


43.99 


41.8 + 3.3 


0.96 + 0.13 


2.52 + 0.07 


028.861+006 


16.91 + 0.94 


103.90 


38.6 + 3.1 


0.36 + 0.06 


1.30 + 0.09 


027.187-008 


15.07 + 0.87 


25.65 


29.3 + 3.1 


0.34 + 0.27 


1.76 + 0.10 


029.436-017 


3.17 + 0.36 


86.19 


36.7+11.0 


0.29 + 0.19 


0.43 + 0.27 


019.923-025 


7.54 + 0.48 


64.30 


46.3 + 7.8 


0.39 + 0.10 


0.89 + 0.16 


024.730+015 


5.65 + 0.44 


109.61 


32.6 + 6.1 


0.25 + 0.16 


0.73 + 0.17 


034.712-059 


9.03 + 0.59 


44.17 


46.4 + 5.5 


0.26 + 0.05 


0.61 + 0.14 


027.925+020 


12.45 + 0.76 


42.65 


34.8 + 2.7 


0.15 + 0.03 


0.52 + 0.11 


037.555+019 


13.67 + 0.83 


86.21 


44.0 + 4.5 


0.47 + 0.08 


1.34 + 0.11 


059.497-023 


17.03 + 0.97 


27.13 


44.4 + 5.1 


0.28 + 0.08 


0.98 + 0.21 


033.384+000 


11.14 + 0.67 


103.88 


40.5 + 4.2 


0.46 + 0.09 


1.29 + 0.10 


025.803-015 


16.41 + 0.94 


92.06 


39.4 + 3.4 


0.36 + 0.06 


1.26 + 0.09 


025.411+010 


12.10 + 0.74 


95.91 


37.7 + 3.2 


0.33 + 0.06 


1.10 + 0.09 


059.355-020 


10.52 + 0.69 


29.53 


49.0 + 5.0 


0.40 + 0.09 


1.04 + 0.17 


025.393+004 


11.27 + 0.66 


-13.02 


34.9 + 2.8 


0.06 + 0.02 


-0.34 + 0.23 


023.385+018 


9.96 + 0.71 


75.75 


37.9 + 4.0 


0.29 + 0.06 


0.89 + 0.11 


Mean 


13.63 






0.35 


1.53 


Standard Deviation 


9.51 






0.19 


1.05 


Median 


11.69 






0.34 


1.01 
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Table 3. HCO+ Properties 



Source 




log(M iso ) 


log(M„ ir ) 




[km s~ x ] 


[M ] 






[M ] 


019.884-053 


5.8 + 0.1 


3.99 + 


.04 


3 


.69 + 0.06 


028.861+006 


5.7 + 0.1 


3.05 + 


.05 


3 


.25 + 0.07 


027.187-008 


7.4 + 0.3 


3.63 + 


.05 


3 


.44 + 0.35 


029.436-017 


3.4 + 0.2 


2.28 + 


.08 


2 


.70 + 0.30 


019.923-025 


4.8 + 0.2 


2.71 + 


.05 


3 


.13 + 0.12 


024.730+015 


3.7 + 0.2 


2.61 + 


.07 


2 


.72 + 0.27 


034.712-059 


9.0 + 0.7 


2.20 + 


.09 


3 


.50 + 0.11 


027.925+020 


4.9 + 0.2 


2.06 + 


.10 


2 


.73 + 0.10 


037.555+019 


8.9 + 0.3 


2.98 + 


.07 


3 


.75 + 0.08 


059.497-023 


3.4 + 0.1 


2.10 + 


.20 


2 


.69 + 0.12 


033.384+000 


4.8 + 0.2 


2.99 + 


.06 


3 


.20 + 0.09 


025.803-015 


3.7 + 0.2 


2.67 + 


.06 


2 


.86 + 0.09 


025.411+010 


6.9 + 0.2 


2.52 + 


.07 


3 


.38 + 0.08 


059.355-020 


2.8 + 0.1 


1.82 + 


.26 


2 


.67 + 0.10 


025.393+004 


7.3 + 0.2 


1.37 + 


.22 


2 


.68 + 0.14 


023.385+018 


3.8 + 0.1 


2.58 + 


.06 


2 


.80 + 0.10 


Mean 




3.09 






3.24 


Standard Deviation 




0.87 






0.41 


Median 




2.60 






3.01 
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Table 4. Photometry 



Source 



A 


Flux 


Beam 


Ref 


H 


[Jy] 


["] 




o 

8 


n a 1 in no 

0.61 ± 0.03 


1 O 

18 


1 


14 


n to i n n/ - * 

0.73 ± 0.06 


18 


1 


21 


o /i o i n nn 

3.43 ± 0.22 


18 


1 


25 


11.07 


300x45 


2 


60 


OA/1 4 r 1 O f-r 

394.45 ± 3.75 


165x87 


3 


65 


oor~ /in i /ii on 

285.49 ± 41.20 


600x60 


4 


70 


a f\ t-7 on i o nn 

407.22 ± 2.99 


62.4 


3 


90 


oil rn i oo a a 

211.59 ± 00.44 


/?nn nn 

600x90 


4 


100 


1155.65 ± 6.07 


195x150 


3 


140 


mo n/ - * i o t op" 

598.96 ± 37.25 


«aa i nn 

700x100 


4 


1 £?n 

160 


Tn c i iic Tn 

795 ± 115.70 


700x150 


A 

4 


OCA 

850 


O O A 


1 o c 

13.5 


5 


i i nn 

1100 


/i n 1 i n o t 

4.91 ± 0.37 


o o 

33 


6 


8 


o o o i n 1 o 

3.28 ± 0.13 


18 


1 


12 


T 1 O I <~i O Z" 4 

7.12 ± 0.36 


18 


1 


14 


15.34 ± 0.94 


18 


1 


21 


34.01 ± 2.04 


18 


1 


o c 

25 


1 f\r on 

105.20 


onn a tz 

300x45 


2 


60 


1945.55 ± 7.08 


165x87 


3 


65 


1/tno n^i i it/" 1 or 

1493.04 ± 176.25 


600x60 


4 


70 


1 r r r il T A A I O TO 

1567.44 ± 3.73 


62.4 


3 


nn 

90 


1207.87 ± 88.65 


oc\c\ nn 

600x90 


4 


100 


oo/^o no i 1 n no 

3263.03 ± 19.92 


195x150 


3 


140 


839.89 ± 2.29 


700x100 


4 


1100 


4.39 ±0.29 


33 


6 


8 


1.35 ±0.06 


18 


1 


12 


3.06 ±0.16 


18 


1 


14 


5.22 ±0.32 


18 


1 


21 


16.77 ± 1.01 


18 


1 


25 


39.20 


300x45 


2 


60 


797.61 ±6.92 


165x87 


3 


65 


527.15 ±98.70 


600x60 


4 


70 


640.54 ±3.38 


62.4 


3 


90 


360.93 ± 120.56 


600x90 


4 


100 


1634.41 ± 14.90 


195x150 


3 


140 


447.09 ± 32.62 


700x100 


4 



019.884-053 



028.861+006 



027.187-008 



-25 - 



Table 4 — Continued 



Source A Flux Beam Ref 1 
[Mm] [Jy] ["] 

160 562.54 ±84.36 700x150 4 

1100 2.20 ±0.22 33 6 

029.436-017 8 1.33 ±0.05 18 1 

12 1.60 ±0.10 18 1 

14 2.33 ±0.15 18 1 

21 3.18 ±0.20 18 1 

65 22.21 ±3.15 600x60 4 

70 18.39 ±4.16 62.4 3 

90 24.28 ±0.65 600x90 4 

160 130.82 ± 22.81 700x150 4 

1100 0.67 ± 0.14 33 6 

019.923-025 8 0.95 ± 0.04 18 1 

12 2.33 ±0.13 18 1 

14 4.15 ±0.25 18 1 

21 5.50 ±0.34 18 1 

25 6.74 300x45 2 

60 13.94 ± 2.27 165x87 3 

70 35.40 ±3.19 62.4 3 

90 32.14 ± 1.06 600x90 4 

100 315.19 ± 23.91 195x150 3 

140 96.56 ±9.75 700x100 4 

160 140.93 ± 22.97 700x150 4 

1100 3.02 ± 0.22 33 6 

024.730+015 8 1.17 ±0.05 18 1 

12 1.52 ±0.08 18 1 

14 1.74 ±0.11 18 1 

21 5.44 ±0.33 18 1 

25 15.20 300x45 2 

60 325.50 90x300 2 

65 177.78 ± 7.83 600x60 4 

70 188.39 ±4.32 62.4 3 

100 1279 180x300 2 

140 649.83 ±62.51 700x100 4 

160 749.97 ± 72.18 700x150 4 

1100 1.39 ± 0.19 33 6 

034.712-059 8 8.13 ±0.33 18 1 
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Table 4 — Continued 



Source A 
H 



Flux 

[Jy] 



10.26 ±0.51 
13.00 ±0.79 
26.89 ± 1.61 

43.70 
290.98 ±2.48 
295.41 ±63.67 

250.61 ±3.20 

178.62 ±65.62 
1086.15 ±9.72 
265.17 ±33.81 

2.46 
0.94 ±0.04 
1.14 ±0.08 
1.68 ±0.11 
6.49 ±0.40 
402.72 ±2.33 
264.64 ±22.53 

218.93 ±3.17 
275.78 ±2.76 

1137.67 ± 10.42 

304.94 ±53.19 
404.62 ±31.68 

1.86 ±0.17 
8.57 ±0.35 
14.83 ±0.74 
24.74 ± 1.51 
37.87 ±2.27 

85.25 
538.36 ±5.57 
380.24 ± 128.67 

510.17±3.26 
288.24 ± 104.53 
937.88 ± 10.59 
429.70 ±86.90 
537.72 ± 143.92 
5.45 



Beam Rcf 1 

["] 



18 


1 


18 


1 


18 


1 


300x45 


2 


165x87 


3 


600x60 


4 


62.4 


3 


600x90 


4 


195x150 


3 


700x100 


4 


33 


6 


18 


1 


18 


1 


18 


1 


18 


1 


165x87 


3 


600x60 


4 


62.4 


3 


600x90 


4 


195x150 


3 


700x1 00 

i UUA1UU 


4 


700x150 


4 


33 


6 


18 


1 


18 


1 


18 


1 


18 


1 


300x45 


2 


165x87 


3 


600x60 


4 


62.4 


3 


600x90 


4 


195x150 


3 


700x100 


4 


700x150 


4 


13.5 


5 



12 
14 
21 
25 
60 
65 
70 
90 
100 
140 
1100 

027.925+020 8 
12 
14 
21 
60 
65 
70 
90 
100 
140 
160 
1100 

037.555+019 8 
12 
14 
21 
25 
60 
65 
70 
90 
100 
140 
160 
850 



Table 4 — Continued 



Source A 
H 



Flux 

[Jy] 



3.11 ±0.25 
0.08 ±0.01 
1.07 ±0.52 

1.65 ±0.11 

4.66 ±0.29 
10.71 

122.40 ± 1.07 
64.53 ±9.10 
33.92 ±3.98 
37.47 ±7.78 
462.16 ±4.80 
96.32 ±39.96 
192.87 ±20.75 
1.37 ±0.23 
0.07 ±0.01 
0.21 ±0.02 
1.22 ±0.07 
5.27 ±0.32 
122.77 ±35.07 

124.36 ±4.99 
145.52 ± 15.97 

279.37 ±76.98 
399.49 ± 143.59 

4.13 ±0.28 
2.26 ±0.09 
8.55 ±0.43 
18.27 ± 1.11 
68.66 ±4.12 
790.50 ±6.81 
349.18 ±35.48 
417.27 ±4.03 

345.45 ± 60.75 
1846.48 ± 16.71 

351.46 ± 120.18 
282.21 ± 135.45 

4.95 



Beam Ref 1 

["] 



33 


6 


18 


1 


18 


1 


18 


1 


18 


1 


300x45 


2 


165x87 


3 


600x60 


4 


62.4 


3 


600x90 


4 


195x150 


3 


700x100 


4 


700x150 


4 


33 


6 


18 


1 


18 


1 


18 


1 


18 


1 


600x60 


4 


62.4 


3 


600x90 


4 


700x100 


4 


700x150 


4 


33 


6 


18 


1 


18 


1 


18 


1 


18 


1 


165x87 


3 


600x60 


4 


62.4 


3 


600x90 


4 


195x150 


3 


700x100 


4 


700x150 


4 


13.5 


5 



1100 

059.497-023 8 
12 
14 
21 
25 
60 
65 
70 
90 
100 
140 
160 
1100 

033.384+000 8 
12 
14 
21 
65 
70 
90 
140 
160 
1100 

025.803-015 8 
12 
14 
21 
60 
65 
70 
90 
100 
140 
160 
850 
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Table 4 — Continued 



Source A 
H 



Flux 

[Jy] 



1.90 ±0.19 
2.33 ±0.10 
3.46 ±0.17 
4.70 ±0.29 
8.43 ±0.51 
17.50 
109.90 
152.96 ± 10.85 
142.10 ±3.36 
75.81 ± 11.84 
157.19 ±43.69 
279.51 ±30.24 
1.01 ±0.16 
1.82 ±0.07 
5.24 ±0.27 
8.72 ±0.53 
25.60 ± 1.54 

41.16 
333.15 ± 1.90 
87.50 ±3.36 
879.54 ±7.21 
187.37 ±40.37 
0.38 ±0.20 
2.29 ±0.09 
4.78 ±0.24 
7.14 ±0.44 
17.00 ± 1.02 
1116.91 ±8.23 
430.67 ± 115.07 
804.15 ±4.61 
307.23 ±59.00 
2350.42 ±26.44 
569.92 ±45.48 
663.00 ±90.06 
2.31 ±0.22 
23.47 ±0.96 



Beam Rcf 1 

["] 



33 


6 


18 


1 


18 


1 


18 


1 


18 


1 


300x45 


2 


90x300 


2 


600x60 


4 


62.4 


3 


600x90 


4 


700x100 


4 


700x150 


4 


33 


6 


18 


1 


18 


1 


18 


1 


18 


1 


300x45 


2 


165x87 


3 


62.4 


3 


1 95x1 50 


O 


700x100 


4 


33 


6 


18 


1 


18 


1 


18 


1 


18 


1 


165x87 


3 


600x60 


4 


62.4 


3 


600x90 


4 


195x150 


3 


700x100 


4 


700x150 


4 


33 


6 


18 


1 



1100 

025.411+010 8 
12 
14 
21 
25 
60 
65 
70 
90 
140 
160 
1100 

059.355-020 8 
12 
14 
21 
25 
60 
70 
100 
140 
1100 

025.393+004 8 
12 
14 
21 
60 
65 
70 
90 
100 
140 
160 
1100 

023.385+018 8 
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Table 4 — Continued 



Source 


A 


r lux 


Beam 


n n f 1 

Kei 




r l 

[Mm] 


r t i 

[Jy] 


[7/1 
[ ] 






12 


39 00 ± 1 95 


18 


1 




14 


50 99 ± 3 11 

'J\J • tj tj _1_ KJ . -L _L 


18 


1 




21 


54 71 ± 3 28 


18 


1 




25 


76.45 


300x45 


2 




60 


281.63 ± 2.84 


165x87 

1 \JKJj\-KJ 1 


3 




65 


262.22 ± 38.11 


600x60 


4 




70 


239.14 ± 3.44 


62.4 


3 




100 


157 45 ± 16 07 

17 1 ■ 1 KJ — 1 — 1 VI . V 1 J 


195x150 


3 




160 


331 91 ± 25 58 

KJ KJ -i- ■ kJ J- — 1 — — » KJ ■ KJ L ' 


700x150 

1 V / V / ^A. 1 KJ\J 


4 




1100 


2.26 ± 0.21 


33 


6 


1 91 30+0 fifi 


12 


i oxnon 


300x45 


2 




25 


21 ± 1.00 


300x45 


2 




60 


•357 _|_ 21 nn 

KJKJ 1 — 1— ± .XJXJ 


90x300 


2 




90 


70ft 59 _|_ OO QQ 
1 UXJ.kjkj —1— oo.uo 


600x90 


4 




100 


685 ± 55 00 


1 80x300 


2 




140 


843 33 ± 70 25 


700x100 


4 




143 


1615 ± 323 00 


105 


7 




185 


2317 ± 463 00 


102 


7 




350 


410 ± 82.00 


120 


8 




450 


40 + no 


18 


g 




800 


« 2 + n 02 


16 


g 




850 


17 ± 3.40 


18 


10 




1100 


5.6 ± 1.10 


18 




Izo.U / -D.ol 


ok 
zo 


1011 nn 
lo ± l.UU 


OUUX40 



z 




60 


330 ± 46.00 


90x300 


2 




65 


263.84 ±73.24 


600x60 


4 




100 


1166 ± 117.00 


180x300 


2 




140 


796.87 ±277.15 


700x100 


4 




160 


784.67 ±319.67 


700x150 


4 




350 


290 ± 58.00 


120 


8 


W3(OH) 


12 


40.58 ±4.06 


300x45 


2 




20 


270 ± 27.00 


49 


8 




25 


670 ± 67.00 


49 


8 




30 


1400 ± 140.00 


49 


8 




40 


4000 ± 400.00 


49 


8 




60 


7000 ± 700.00 


49 


8 
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Table 4 — Continued 



Source A 
[Mm] 



65 
70 
100 
140 
160 
180 
200 
250 
300 
350 
500 
600 
800 
1000 
S231 12 
25 
60 
100 
160 
350 

S252A 12 
25 
60 
100 
140 
350 

RCW142 12 
25 
60 
65 
90 
100 
160 
350 

W28A2 12 
25 



Flux 

[Jy] 



6918.19 ± 198.06 
8500 ±850.00 
9000 ± 900.00 
6700 ±670.00 
8645.55 ± 1006.47 
4900 ±490.00 
4100 ±410.00 
2400 ± 240.00 
1400 ± 140.00 
1130 ±230.00 
250 ±25.00 
135 ± 14.00 
51 ±5.10 
24 ±2.40 
5.6 ±0.20 
75 ±3.70 
722 ±72.00 
1310 ± 131.00 
1841.75 ±94.64 
522 ± 100.00 
16 ±0.60 
77 ±3.00 
10321 ±34.00 
1715 ± 189.00 
844.94 ± 131.43 
320 ±64.00 

< 42 

< 281 
5476 ±986.00 

2623.98 ±59.43 
1616.73 ± 134.92 
13129 ± 1313.00 
8480.46 ± 824.94 

670 ± 130.00 
199 ± 12.00 

2190 ± 131.00 



Beam Ref 1 
["] 



600x60 


4 


49 


8 


49 


8 


49 


8 


700x150 


4 


49 


8 


49 


8 


49 


8 


49 


8 


120 


8 


49 


8 


49 


8 


49 


8 


49 


8 


300x45 


2 


300x45 


2 


90x300 


2 


180x300 


2 


700x150 


4 


120 


8 


300x45 


2 


300x45 


2 


90x300 


2 


180x300 


2 


700x100 


4 


120 


8 


300x45 


2 


300x45 


2 


90x300 


2 


600x60 


4 


600x90 


4 


180x300 


2 


700x150 


4 


120 


8 


300x45 


2 


300x45 


2 
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Table 4 — Continued 



Source 



A 


Flux 


Beam 


Ref 


H 


[Jy] 


["] 




fin 


1 97Q0 -i- ^1 ne. nn 
iz/yu it oiyo.uu 


onv^nn 
yuxouu 


9 
Z 


1UU 


o^7Qn J- fifiorc nn 
ZO/oU ± ooyo.uu 


loUxoUU 


Z 


Q^n 

oOU 


1 rccn -u Qon nn 

lOoU zt oZU.UU 


1 on 
IzU 


Q 
O 


1 


llo.O zt 11. ou 


OUUX40 


9 

z 


0£ 


oqq 4- 17 nn 
Zoy ± 1/ .uu 


OUUX40 


z 


fin 
ou 


1 fii 1 _i_ oofi nn 
1D1 1 it ZZO.UU 


on-^^nn 
yuxouu 


z 


fi£ 
DO 


DOo.41 zt 01U.U1 


ouuxou 


/] 

4 


fiQ 

oy 


ZOUU 


KA 
04 


1 
1Z 


on 
yu 


cn^ -4- Qc;r; qi 
OUO.ZO zt oOO.ol 


ouuxyu 


4 


1 nn 

1UU 


0700 1 (ZQa nn 
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1 REFERENCES : (1) MSX Database. (2) IRAS PSC 1988. (3) Mottram et al. 2010. (4) AKARI Database. 
(5) SCUBA Legacy Catalogs. (6) Aguirre et al. 2011. (7) Mookerjea et al. 1999. (8) Mueller ct al. 2002. 
(9) Dent, Matthews, & Ward-Thompson 1998. (10) Jenness, Scott, & Padman 1995. (11) McCutcheon et 
al. 1995. (12) Chini, Henning, & Pfau 1991. (13) Tothill 1999. (14) Chini et al. 1986a. (15) Fazio et al. 
1978. (16) Chini et al. 1986b. (17) Jaffe et al. 1984. (18) Jenness 1996. (19) Evans et al. 1979. (20) Dyck 
& Simon 1977. (21) Sticr ct al. 1984. (22) Gucrtlcr ct al. 1991. (23) Cheung et al. 1980. (24) Soifer & 
Pipher 1975. (25) Hatchell et al. 2000. (26) Zhou et al. 1994. (27) Ellis et al. 1990. (28) Evans ct al. 
1981b. (29) Moriarty-Schieven, Snell, & Hughes 1991. (30) Gordon 1990. (31) Werner et al. 1979. 



Table 5. Bayesian linear correlations 
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Fig. 1. — HCO + (3-2) maps for BGPS/RMS sources. The lowest intensity contour and 
the contour interval are both 2 a in each map. The percentage of the 2 a intensity is also 
specified. 
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Fig. 2. — Histograms of source properties. The shaded histograms represent the BGPS/RMS 
sources and the bold line histograms represent the water-maser sources. 
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Fig. 3.- (a) log(L')-log(M iso ) for HCO+ (3-2). (b) log(L')-log(M Mr ) for HCO+ (3-2). (c) 
log(L 6o ;)-log(L') for HCO+ (3-2). (d) log(L 6oi )-log(L') for HCN (1-0). The inset histograms 
display the range of slopes determined from the Bayesian routine. 



